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Abstract 



We theoretically study electronic states in graded-gap junctions of IV- VI 
compounds with band inversion. Using a two-band model within the k ■ p 
approximation and assuming that the gap and the gap centre present linear 
profiles, we demonstrate the existence of a set of localized states along the 
growth direction with a discrete energy spectrum. The envelope functions 
are found to be combination of harmonic oscillator eigenfunctions, and the 
corresponding energy levels are proportional to the square root of the quantum 
number. The level spacing can be directly controlled by varying the structure 
thickness. 
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Narrow-gap IV- VI compounds like Pbi-zSn-jTe and Pbi-^Sn^Se present band inversion 
under compositional variation. In a band-inverted heteroj unction the fundamental gap has 
opposite signs on each side [[!]]. For instance, Pbi_ x Sn x Se undergoes band inversion as the 
Sn molar fraction is increased: At x — 0.14 the gap vanishes whereas the negative one 
reaches the magnitude of the PbSe gap at x = 0.28. Recently, such heteroj unctions have 
received much attention because subbands of electron-like and hole-like localized interface 
states are formed with energy lying within the fundamental gap f2)-|5|], a particular feature 
not met in most common III-V heterojunctions. In those compounds the electronic states 
near the gap are properly described by means of a two-band model using the effective k ■ p 
approximation 0, where quite strong coupling of host bands in the semiconductor takes 
place. The equation governing conduction- and valence-band envelope-functions in a simple 
two-band model, neglecting far-band corrections, is a Dirac-like equation. In view of the 
analogy existing between the two-band model and the Dirac equation, the exact solution 
can be found since one can use elaborated techniques like those related to supersymmetric 
quantum mechanics |J. 

In this paper we exploit further such a formal similarity to demonstrate for the first time 
the existence of a new type of localized electronic states whose envelope functions resemble 
those of the harmonic oscillator but, unlike the case of the well-known parabolic quantum 
wells, there are no classical turning points. In the present structure the confining poten- 
tial leading to localized states enters in the equation of motion as a position-dependent 
mass term, whereas in the case of parabolic quantum wells the conduction-band modula- 
tion appears as an electrostatic-like interaction within the framework of a one-band model 
Hamiltonian. In our approach, we require a graded-gap structure with band inversion. An 
appropriate graded doping may create a modulation of both conduction- and valence-band 
with linear profiles, which manifest themselves through the occurrence of a linear scalar 
potential in the Dirac-like Hamiltonian. For our present purposes, we then take advan- 
tage that there exists a number of linear potentials for which the Dirac equation reduces to 
Schrodinger-like harmonic oscillator equations f7j, and then exact solutions can be found in 
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a closed form. 

The two-band model Hamiltonian in the absence of external fields is of the form 



where the z axis is perpendicular to the hetero junction (it is assumed that the growth 
direction is [111]), E g (z) stands for the position dependent gap, a y , a z and (3 are the usual 
4x4 Dirac matrices, v± and v z have dimensions of velocity and they are related to the 
Kane's matrix elements. As usual, it is assumed that these matrix elements are constant 
through the whole heterostructure. Since the gap and the gap centre depend only upon 
z, the transversal momentum is a constant of motion and we can set the y axis parallel 
to this component. In a graded-gap structure of thickness L with band inversion we have 
E g (z) = Kz, where K = (E gR — E 9 l)/L. Here E gR > [E 9 l < 0) is the magnitude of the 
gap at z R (—z l ) with L = z R + z L and E gR /E gL = -z R /z L . 

In the two-band model there are four envelope-functions including spin, and we arrange 
them in a four component spinor F(r), whose upper and lower components give the coeffi- 
cients of the Lq and Lg" parts of the wave function. This spinor satisfies the equation 



where V(z) gives the position of the gap centre. The way Viz) changes from one material to 
another is not well understood, so that it is often considered that the misalignment follows 
the same profile of E g {z) §. Therefore V(z) = [\/2)E g {z) where A/2 = (V R - V L )/{E gR - 
E 9 l), V r and Vl being the gap centre at z R and —zl, respectively. We shall see below 
that the existence of localized states requires |A| < 1, namely the gaps must overlap (type I 
hetero junctions). The same condition is found for the existence of interface states in band 
inverted junctions ||. As we have already mentioned, the momentum perpendicular to the 
interface is conserved, and therefore we look for solutions of the form F(r) = F(z) exp(i r_|_ • 
p±/h) to Eq. (0). The function F(z) satisfies the following equation 



U = a y v±p ± + a z v z p z + - PE g (z), 



(1) 



HF(r) = [E-V(z)]F(r], 



(2) 
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A simple way to solve this equation is the Feynman-Gell-Mann ansatz 



F(z) 



a y v±p± + a z v z p z + - /3E g (z) + E — V(z) ) x( z )- 



(4) 



After a little algebra we have 
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dz 2 % 2 v 2 



E g {zY-[E -V{z)Y+ vipi 



+ M\x(z) = 0, 



(5) 



where the 4x4 constant matrix M is given by M = —i(K/2hv z )a z (j3 — \), whose eigenvalues 
±/i = ±{K/2%v z )y/l - A 2 are real since we take |A| < 1. Setting x( z ) = f±{ z )4 > ±i where <f>± 
are the eigenvectors corresponding to the eigenvalues ±/i, we obtain 

E 2 - v\p\ 



d 2 2 2 \KE 
— + fx z + -prr z 
dz 2 h v 2 



h 2 v 2 



=F/i 



f±( z ) 



(6) 



which clearly reduces to a nonrelativistic oscillator equation by carrying out a suitable 
translation of the origin of coordinates (recall that fi 2 > 0). Thus, inverting the various 
transforms needed to arrive at ([]), it is not difficult to demonstrate that envelope functions 
are simply combinations of Hermite polynomials times a decreasing exponential factor similar 
to that of the harmonic oscillator, provided that f(z) vanishes at \z\ — > oo in a suitable way. 
The corresponding bound levels can also be found in a simple fashion, and they are given 
by 



Ei 



nhv z KVl — A 2 + v \p\ 



(7) 



n being a positive integer. Notice that energy levels increase as the square root of K and n, 
thus being no longer equally spaced. This behaviour is usually obtained in "relativistic-like" 
oscillator equations 0. As an example, let us consider the case of selenides with symmetric 
band inversion, i.e., — E 9 l = E g R = E g and Vl = Vr (A = 0). Typical parameters are 
Eg = 0.15 eV, hv z = 2.76 eV. For L = 2zr = 2zl = 600 A and p± = we obtain four bound 
state levels with energy less than E g /2: Ex = 37.2 meV, E 2 = 52.5 meV, E 3 = 64.3 meV 
and £4 = 74.3 meV. Hence the level spacings are E2 — E\ = 15.3 meV, £3 — E2 = 11.8 meV, 
£4 — E 3 = 10.0 meV, showing that levels are not equally spaced but deviation from that 
behaviour is actually small. 



To summarize, we have discussed in some detail the electronic structure of graded-gap 
hetrojunctions of IV- VI compuounds presenting band inversion. The theoretical analysis 
is based on the two-band model arising in the k • p approach which, neglecting far-band 
couplings, becomes completely analogous to a Dirac-like equation with linearly rising scalar- 
and electrostatic-like potentials. The Dirac Hamiltonian for linear potentials is exactly 
solvable, leading to envelope functions that can be expressed as a sum of harmonic oscillator 
functions. Therefore, carriers are spatially localized close to z = 0, i.e., the plane where 
band inversion occurs. However, notice that the maximum value of envelope functions is 
located at zq = —2XE jK\pY — A 2 since Eq. (|j) is a harmonic oscillator equation centered 
at Zq. This is to be compared to interface states lying within the fundamental gap which, 
in absence of external fields, are centered at the crossing point of Lg and Lg" bands [[|. In 
addition, we have found that energy levels increase as the square root of effective coupling 
K and the quantum number n. Let us stress that the level spacing and spatial extend of 
envelope functions depend on the value of K, and as a consequence they can be controlled 
by varying the values of E g i and E g R as well as the thickness of the structure L. Finally, 
some words concerning the validity of the present model are in order. It is known that 
the Hamiltonian (|I|) has limitations since it neglects far-band corrections. Those effects can 
be evaluated by means of the standard second-order perturbation theory ||, although we 
do not attempt to carry out such computation here. However, we can confidently expect 
that our results are valid, at least qualitatively, since in most selenides such corrections 
only cause small deviations of the results predicted from ([J). In fact, interface states lying 
within the fundamental gap remain even if far-band terms are included in the Hamiltonian, 
while only minor modifications of the dispersion relation are observed ||. We hope that our 
results may encourage experimental effort in this field for two reasond. First, to validate or 
discard the existence of localized states in actual band-inverted junctions other than interface 
states. Second, the feature of having oscillator-like states along with their associate discrete 
spectrum may be the basis for designing new devices and applications. 
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